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We consider a Bose-Einstein condensate of ultracold atoms loaded into a square optical lattice 
and subject to a static force. For vanishing atom-atom interactions the atoms perform periodic 
Bloch oscillations for arbitrary direction of the force. We study the stability of these oscillations 
for non- vanishing interactions, which is shown to depend on an alignment of the force vector with 
respect to the lattice crystallographic axes. If the force is aligned along any of the axes, the mean 
field approach can be used to identify the stability conditions. On the contrary, for a misaligned 
force one has to employ the microscopic approach, which predicts periodic modulation of Bloch 
oscillations in the limit of a large forcing. 



A Bose-Einstein condensate (BEC) in optical lattices 
has intrigued a rapidly growing interest as it provides 
an experimentally realizable system with controllable in- 
teractions. A variety of phenomena concerning differ- 
ent aspects of physics for quantum many-body systems 
has been studied, such as the superfluid-Mott insulator 
quantum phase transition 0, Q , BEC-BCS crossover for 
fermionic gases Q , and quantum transport in accelerated 
lattices, where atoms exhibit fundamental quantum ef- 
fects such as the Wannier-Stark ladder 4J, Landau-Zener 
tunnelling [5[ , and Bloch oscillations [J, ||| - phenomena 
usually associated with electron in solid crystall. This 
work deals with the last mentioned problem, namely, 
Bloch oscillations (BO) of condensed atoms in optical 
lattices. We would like to mention that besides pure aca- 
demic interest this problem also has an applied aspect 
because BO provide a tool for precision measurement of 
gravitational field and inter- atomic interaction constant. 

Until quite recently almost all theoretical, numerical 
and experimental studies of BO concerned ID or quasi 
ID lattices (see Ref. @, IE f° r the contemporary re- 
views). Nowadays one observes a growing interest in BO 
in multidimensional lattices [13, [ll, 12, 13 , 14, 1^|. In the 
single-particle approach this problem was considered in 
Refs. [10|,lll|. It was shown that an increase of the lattice 
dimensionality introduces new effects not present in the 
ID lattice. Some predictions of these works were later 
on confirmed in the experiment with the array of opti- 
cal guides [13| , where one uses a formal analogy between 
the Maxwell and Schrodinger equations. The experiment 
with a BEC of interacting atoms addresses the further 
questions [lij], in particular, the question about the sta- 
bility of multidimensional BO. Indeed, it is known that 
a BEC in optical lattices can be dynamically unstable, 
which quantum-mechanically means decoherence of the 
BEC [161 ] . In the present work we study the conditions 
under which the dynamical instability is suppressed and, 
hence, multidimensional BO are stable. Unlike ID lat- 
tices, these conditions are shown to involve an alignment 



of the static force vector with respect to the crystallo- 
graphic axes of the lattice. We also argue in the work 
that by changing the angle between the primary lattice 
vectors and the force vector one may observe a transition 
from the mean-field to the microscopic Bloch dynamics. 

To simplify the equations we shall consider the two- 
dimensional case throughout the paper, - generalization 
of the results in three dimensions is straightforward. The 
Bose-Hubbard Hamiltonian of atoms in the tilted 2D lat- 
tice reads, 
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where J x%v are the hopping matrix elements in x and 
y directions, W microscopic atom-atom interaction con- 
stant, d lattice period, and F X y the projections of the 
static force vector on the lattice axes. The Hilbert space 
of (jTJ) is spanned by the Fock states |n) = \n m j), where 
Sm i n ™.,i = N - the total number of atoms. Since in the 
coordinate representation the Fock states are given by 
the symmetrized product of the localized Wannier func- 
tions, we shall refer to this basis as the Wannier basis. 
The translational invariance of the system, broken by the 
static term, can be actually recovered by using the gauge 
transformation [2J. Then the Hamiltonian (fl]) takes the 
form 
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where 



= dF x .y/h are the Bloch frequencies associ- 



ated with x and y component of the static force. We also 
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note that in stead of the Wannier basis one can use the 
quasimomentum Fock basis |q) = \q p ,k) for the Hamil- 
tonian @, which we shall refer to as the Bloch basis. 
(Needless to say that in the coordinate representation 
the quasimomentum Fock states are given by the sym- 
metrized product of the extended Bloch functions.) For- 
mally this corresponds to the canonical transformation 
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which implicitly assumes the periodic boundary condi- 
tions. 




Misaligned Force — We begin with the case of a strong 
misaligned force dF x ,dF y 3> J x ,y > W. In order to il- 
luminate situation, we model BO in a small 2D lattice 
by numerically solving the time-dependent Schrodinger 
equation with the Hamiltonian (|2|) for the specified ini- 
tial conditions. As those we consider the superfluid state 
with f/p.fc = NS P fi5k,Q , which approximates the ground 
state of the system for F = and W < J x . y . (Substi- 
tution of this state by the exact ground state practically 
does not affect the final result.) Fig. Q] shows the nu- 
merical results for a 3 x 3 lattice with 7 atoms inside. 
The lower panel in Fig. [1] depicts the mean energy of 
the system, the upper and middle panels show the order 
parameters e x (t) and e y (t) defined as [HI 
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[By replacing the operators in the bracket in Eq. with 
a m i one obtains a similar expression for the or- 
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der parameter e y {t).\ It is seen in the figure that BO 
persist in time but are modulated with some characteris- 
tic period. We would also like to mention that the Bloch 
dynamics displayed in Fig. [1] is converged in the ther- 
modynamic limit, i.e., for given n = N/L 2 the further 
increase of the system size affects neither the modulation 
period nor the shape of modulation. 

To prove that BO in the misaligned lattice are stable 
in the limit of strong forcing and to identify the mod- 
ulation period we proceed as follows. First we intro- 
duce the new wave function \^(t)) through the relation 
|*(t)) = U (t)\^(t)), where U (t) is the evolution oper- 
ator for vanishing atom- atom interactions. The function 
\ty(t)) obviously obeys the equation, 
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On the other hand, the explicit form of the evolution op- 
erator is given by Uo(t) — T'D(t)T, where the unitary 
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FIG. 1: Bloch oscillations of condensed atoms in the 2D lat- 
tice: the order parameter e x (t) (top), e y (t) (middle), and 
the mean energy (bottom). The system parameters are 
N = 7, L = 3 (periodic boundary conditions), J x — J y — J, 
W = 0.2 J, Fd = 20 J, and F/F = {</2/E, The time 

is measured in units of the tunnelling period Tj — 2nh/J. 



operator T represents the transformation from the Wan- 
nier basis |n) to the Bloch basis |q) and the matrix of 
the operator D(t) is diagonal in the Bloch basis, 
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Note that the operator ([5]) tends to the identity operator 
for F x ,F y — > oo. Substituting Uo(t) in Eq. (|4]) by iden- 
tity matrix and noting that the interaction energy opera- 
tor is diagonal in the Wannier basis with integer entries, 
( n l E m ,/ n m ,i{n m ,i - l)|n) = Y, m ,l n L,i~ N ^ ^e conclude 
that the time evolution of the wave function is pe- 

riodic with the period Tw = 2ttH/W [171 ] ■ Coming back 
to the original wave function this result means the peri- 
odic modulation of BO with the frequency low = W/h. 
It is worth stressing that the above proof assumes both 
F x and F y to be large and, hence, the case of aligned 
lattices is excluded. 

Aligned Force — Next we consider the situation where 
the force is aligned along one of the crystallographic axes 
(to be certain, the y-axis in what follows). Within the 
single-particle approach the static force F y would localize 
the atoms in the y-direction. Thus one may expect that if 
F y is large enough the atoms form separate BECs in the 
planes perpendicular to the force vector, weakly coupled 
together as a one-dimensional BEC chain. Introducing 
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new operators Ai ■ 
Hamiltonian reads 

H eff = 



1 
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and Aj, the effective 
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where VFe// = VF/L. Thus we have reduced the 2D prob- 
lem to a ID problem with the renormalized interaction 
constant. (If one considers 3D lattices, the renormaliza- 
tion is W e ff = W/L 2 .) Moreover, since the mean num- 
ber of atoms N in any site of the effective ID system 
is given by nL, the occupation numbers will be macro- 
scopically large in the thermodynamic limit N, L — > oo, 
n = N/L 2 = const, which justifies the mean field ap- 
proach. 

The mean-field Hamiltonian of the system © reads 
(up to the irrelevant constant terms proportional to 
EiNi = N) 



tAt +1 A l + h.c.) + ^\ A 
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where Ai and A* are pairs of the canonically conju- 
gated variables and the macroscopic interaction constant 
g = W e ffN = Wn. Within the mean-field approach 
the border between stable and unstable (decaying) BO is 
know exactly [19t l2Cl|. Namely, for J/Fd > 0.5 the criti- 
cal value of nonlinearity is a linear function of the static 
force magnitude, while for J/Fd < 0.5 it additionally 
depends on the value of the hopping matrix elements: 
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Obviously, for a fixed nonlinearity g the condition (|8|) can 
also be formulated as a condition on the critical magni- 
tude F cr of the static force. 

The microscopic analysis of BO in the aligned lattice 
confirms our working hypothesis. Choosing the parame- 
ters in such a way that the ID mean-field BO are stable, 
we simulate BO of N — 7 atoms in the 2D lattice with 
L = 3. The dashed line in the upper panel of Fig. [2] 
shows the dynamics of the order parameter e x (t). It is 
seen that e x (t) ~ —1, - therefore we indeed have in-plane 
BECs. We also note that the decay and revival of the or- 
der parameter e y (t) in the lower panel is an artifact due 
to the finite size of our lattice. Indeed, it can be shown 
that the time evolution of e y {t), calculated on the basis 
of the effective Hamiltonian ([6]), obeys the equation 21 1 
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FIG. 2: Dynamics of the order parameters for F/F = (0, 1), 
dashed lines, and F/F = (0.001, VI - 0.001), solid lines. The 
other parameters are the same as in Fig. 1. 



Because W e ff = W/L and N = nL, one has e y (t) = — 1 
in the thermodynamic limit. 

The above analysis of BO in the aligned lattice relies on 
the reduction of a two-dimensional system to an effective 
one-dimensional mean-field problem. It should be espe- 
cially stressed that this reduction is possible only if the 
dynamics of the reduced system is stable. If we choose 
the parameters in the unstable regime the situation be- 
comes totally different. Figure [3] shows the numerical 
results for F = 0.2J/d < F cr , where one-dimensional BO 
suffer from dynamical instability. Unlike in the stable 
regime, BO along y direction now excite the transverse 
degree of freedom and we observe decay of the both or- 
der parameters towards zero. Thus no reduction to one 
dimension is possible. 

Slightly Misaligned Force — Finally we briefly analyze 
an experimentally important situation of a small mis- 
match between the lattice axis and the static field vector, 
i.e., F x <§; F y . The solid lines in Fig. O show the order 
parameters for the same dF — 20J but F x — 0.001F. 
Compared to the case F x = (dashed lines in Fig. [2]), we 
observe the destruction of BEC after time t* w 12Tj. 
This critical time can be understood in terms of the 
mean-field approach as well. Indeed, it is known that 
a stationary BEC is unstable for the quasimomentum k 
outside the first quarter of the Brillouine zone. Since the 
static force causes the linear growth of the quasimomen- 
tum, K x , y (t) = K x<y + F Xt yt/h, the system always enter 
the instability region of the Brillouine zone. However, if 
the static force is strong enough, the system passes the 
instability region so quickly that it 'has no time' to de- 
cay. [In fact this is a physical argument behind Eq. ((5J).] 
In the considered example the strong static force ensures 
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(not mean-field) regime of BO, where they are modulated 
with the frequency defined by the microscopic interaction 
constant. 

Fruitful discussions with A. Buchleitner and finan- 
cial support by Deutsche Forshungsgemeinschaft under 
SPP1116 are gratefully acknowledged. 
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FIG. 3: Bloch dynamics for dF = 0.2J and F/F = (0,1). 
The other parameters are the same as in Fig. 1. 



the fast driving along y direction but simultaneously it 
slowly brings the system to k x = n/2d along x direc- 
tion. As soon as this border of instability is reached 
(t*/Tj = J/4dF x ), we observe an irreversible decay of 
the order parameters. 

Conclusion. In summary, we have studied BO of a 
BEC of atoms in a square lattice for both aligned and 
misaligned static forces. It is shown that in the case 
of aligned force the system may be reduced to a one- 
dimensional chain of mini BECs, which we treat by using 
the mean-field approach. Then the stability diagram of 
this effective ID system defines the critical magnitude 
of the static force above which BO are stable, with no 
excitations of the transverse degrees of freedom. On the 
contrary, in the unstable regime, F < F cr , BO induced 
by the static force excite the transverse modes and, as 
a consequence, one observes BEC distruction and decay 
of BO. Our studies also illuminated importance of the 
alignment. The strong (F > F cr ) but slightly misaligned 
force is shown to slowly intrigue the transverse modes, 
which destabilize BO after some well-defined transient 
time. However, if misalignment is large, BO appear to 
be stable again. This case corresponds to the quantum 
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